3.285 \(\int (a+b x^n) (c+d x^n)^3 \, dx\)

Optimal. Leaf size=99 \[ \frac{c^2 x^{n+1} (3 a d+b c)}{n+1}+\frac{d^2 x^{3 n+1} (a d+3 b c)}{3 n+1}+\frac{3 c d x^{2 n+1} (a d+b c)}{2 n+1}+a c^3 x+\frac{b d^3 x^{4 n+1}}{4 n+1} \]

[Out]

a*c^3*x + (c^2*(b*c + 3*a*d)*x^(1 + n))/(1 + n) + (3*c*d*(b*c + a*d)*x^(1 + 2*n))/(1 + 2*n) + (d^2*(3*b*c + a*
d)*x^(1 + 3*n))/(1 + 3*n) + (b*d^3*x^(1 + 4*n))/(1 + 4*n)

________________________________________________________________________________________

Rubi [A]  time = 0.0727884, antiderivative size = 99, normalized size of antiderivative = 1., number of steps used = 2, number of rules used = 1, integrand size = 17, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.059, Rules used = {373} \[ \frac{c^2 x^{n+1} (3 a d+b c)}{n+1}+\frac{d^2 x^{3 n+1} (a d+3 b c)}{3 n+1}+\frac{3 c d x^{2 n+1} (a d+b c)}{2 n+1}+a c^3 x+\frac{b d^3 x^{4 n+1}}{4 n+1} \]

Antiderivative was successfully verified.

[In]

Int[(a + b*x^n)*(c + d*x^n)^3,x]

[Out]

a*c^3*x + (c^2*(b*c + 3*a*d)*x^(1 + n))/(1 + n) + (3*c*d*(b*c + a*d)*x^(1 + 2*n))/(1 + 2*n) + (d^2*(3*b*c + a*
d)*x^(1 + 3*n))/(1 + 3*n) + (b*d^3*x^(1 + 4*n))/(1 + 4*n)

Rule 373

Int[((a_) + (b_.)*(x_)^(n_))^(p_.)*((c_) + (d_.)*(x_)^(n_))^(q_.), x_Symbol] :> Int[ExpandIntegrand[(a + b*x^n
)^p*(c + d*x^n)^q, x], x] /; FreeQ[{a, b, c, d, n}, x] && NeQ[b*c - a*d, 0] && IGtQ[p, 0] && IGtQ[q, 0]

Rubi steps

\begin{align*} \int \left (a+b x^n\right ) \left (c+d x^n\right )^3 \, dx &=\int \left (a c^3+c^2 (b c+3 a d) x^n+3 c d (b c+a d) x^{2 n}+d^2 (3 b c+a d) x^{3 n}+b d^3 x^{4 n}\right ) \, dx\\ &=a c^3 x+\frac{c^2 (b c+3 a d) x^{1+n}}{1+n}+\frac{3 c d (b c+a d) x^{1+2 n}}{1+2 n}+\frac{d^2 (3 b c+a d) x^{1+3 n}}{1+3 n}+\frac{b d^3 x^{1+4 n}}{1+4 n}\\ \end{align*}

Mathematica [A]  time = 0.105617, size = 90, normalized size = 0.91 \[ \frac{b x \left (c+d x^n\right )^4-x \left (\frac{3 c^2 d x^n}{n+1}+c^3+\frac{3 c d^2 x^{2 n}}{2 n+1}+\frac{d^3 x^{3 n}}{3 n+1}\right ) (b c-a d (4 n+1))}{4 d n+d} \]

Antiderivative was successfully verified.

[In]

Integrate[(a + b*x^n)*(c + d*x^n)^3,x]

[Out]

(b*x*(c + d*x^n)^4 - (b*c - a*d*(1 + 4*n))*x*(c^3 + (3*c^2*d*x^n)/(1 + n) + (3*c*d^2*x^(2*n))/(1 + 2*n) + (d^3
*x^(3*n))/(1 + 3*n)))/(d + 4*d*n)

________________________________________________________________________________________

Maple [A]  time = 0.009, size = 104, normalized size = 1.1 \begin{align*} a{c}^{3}x+{\frac{b{d}^{3}x \left ({{\rm e}^{n\ln \left ( x \right ) }} \right ) ^{4}}{1+4\,n}}+{\frac{{c}^{2} \left ( 3\,ad+bc \right ) x{{\rm e}^{n\ln \left ( x \right ) }}}{1+n}}+{\frac{{d}^{2} \left ( ad+3\,bc \right ) x \left ({{\rm e}^{n\ln \left ( x \right ) }} \right ) ^{3}}{1+3\,n}}+3\,{\frac{cd \left ( ad+bc \right ) x \left ({{\rm e}^{n\ln \left ( x \right ) }} \right ) ^{2}}{1+2\,n}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((a+b*x^n)*(c+d*x^n)^3,x)

[Out]

a*c^3*x+b*d^3/(1+4*n)*x*exp(n*ln(x))^4+c^2*(3*a*d+b*c)/(1+n)*x*exp(n*ln(x))+d^2*(a*d+3*b*c)/(1+3*n)*x*exp(n*ln
(x))^3+3*c*d*(a*d+b*c)/(1+2*n)*x*exp(n*ln(x))^2

________________________________________________________________________________________

Maxima [F(-2)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Exception raised: ValueError} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+b*x^n)*(c+d*x^n)^3,x, algorithm="maxima")

[Out]

Exception raised: ValueError

________________________________________________________________________________________

Fricas [B]  time = 1.59923, size = 695, normalized size = 7.02 \begin{align*} \frac{{\left (6 \, b d^{3} n^{3} + 11 \, b d^{3} n^{2} + 6 \, b d^{3} n + b d^{3}\right )} x x^{4 \, n} +{\left (3 \, b c d^{2} + a d^{3} + 8 \,{\left (3 \, b c d^{2} + a d^{3}\right )} n^{3} + 14 \,{\left (3 \, b c d^{2} + a d^{3}\right )} n^{2} + 7 \,{\left (3 \, b c d^{2} + a d^{3}\right )} n\right )} x x^{3 \, n} + 3 \,{\left (b c^{2} d + a c d^{2} + 12 \,{\left (b c^{2} d + a c d^{2}\right )} n^{3} + 19 \,{\left (b c^{2} d + a c d^{2}\right )} n^{2} + 8 \,{\left (b c^{2} d + a c d^{2}\right )} n\right )} x x^{2 \, n} +{\left (b c^{3} + 3 \, a c^{2} d + 24 \,{\left (b c^{3} + 3 \, a c^{2} d\right )} n^{3} + 26 \,{\left (b c^{3} + 3 \, a c^{2} d\right )} n^{2} + 9 \,{\left (b c^{3} + 3 \, a c^{2} d\right )} n\right )} x x^{n} +{\left (24 \, a c^{3} n^{4} + 50 \, a c^{3} n^{3} + 35 \, a c^{3} n^{2} + 10 \, a c^{3} n + a c^{3}\right )} x}{24 \, n^{4} + 50 \, n^{3} + 35 \, n^{2} + 10 \, n + 1} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+b*x^n)*(c+d*x^n)^3,x, algorithm="fricas")

[Out]

((6*b*d^3*n^3 + 11*b*d^3*n^2 + 6*b*d^3*n + b*d^3)*x*x^(4*n) + (3*b*c*d^2 + a*d^3 + 8*(3*b*c*d^2 + a*d^3)*n^3 +
14*(3*b*c*d^2 + a*d^3)*n^2 + 7*(3*b*c*d^2 + a*d^3)*n)*x*x^(3*n) + 3*(b*c^2*d + a*c*d^2 + 12*(b*c^2*d + a*c*d^
2)*n^3 + 19*(b*c^2*d + a*c*d^2)*n^2 + 8*(b*c^2*d + a*c*d^2)*n)*x*x^(2*n) + (b*c^3 + 3*a*c^2*d + 24*(b*c^3 + 3*
a*c^2*d)*n^3 + 26*(b*c^3 + 3*a*c^2*d)*n^2 + 9*(b*c^3 + 3*a*c^2*d)*n)*x*x^n + (24*a*c^3*n^4 + 50*a*c^3*n^3 + 35
*a*c^3*n^2 + 10*a*c^3*n + a*c^3)*x)/(24*n^4 + 50*n^3 + 35*n^2 + 10*n + 1)

________________________________________________________________________________________

Sympy [A]  time = 1.7615, size = 1540, normalized size = 15.56 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+b*x**n)*(c+d*x**n)**3,x)

[Out]

Piecewise((a*c**3*x + 3*a*c**2*d*log(x) - 3*a*c*d**2/x - a*d**3/(2*x**2) + b*c**3*log(x) - 3*b*c**2*d/x - 3*b*
c*d**2/(2*x**2) - b*d**3/(3*x**3), Eq(n, -1)), (a*c**3*x + 6*a*c**2*d*sqrt(x) + 3*a*c*d**2*log(x) - 2*a*d**3/s
qrt(x) + 2*b*c**3*sqrt(x) + 3*b*c**2*d*log(x) - 6*b*c*d**2/sqrt(x) - b*d**3/x, Eq(n, -1/2)), (a*c**3*x + 9*a*c
**2*d*x**(2/3)/2 + 9*a*c*d**2*x**(1/3) + a*d**3*log(x) + 3*b*c**3*x**(2/3)/2 + 9*b*c**2*d*x**(1/3) + 3*b*c*d**
2*log(x) - 3*b*d**3/x**(1/3), Eq(n, -1/3)), (a*c**3*x + 4*a*c**2*d*x**(3/4) + 6*a*c*d**2*sqrt(x) + 4*a*d**3*x*
*(1/4) + 4*b*c**3*x**(3/4)/3 + 6*b*c**2*d*sqrt(x) + 12*b*c*d**2*x**(1/4) + b*d**3*log(x), Eq(n, -1/4)), (24*a*
c**3*n**4*x/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 50*a*c**3*n**3*x/(24*n**4 + 50*n**3 + 35*n**2 + 10*n +
1) + 35*a*c**3*n**2*x/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 10*a*c**3*n*x/(24*n**4 + 50*n**3 + 35*n**2 +
10*n + 1) + a*c**3*x/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 72*a*c**2*d*n**3*x*x**n/(24*n**4 + 50*n**3 + 3
5*n**2 + 10*n + 1) + 78*a*c**2*d*n**2*x*x**n/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 27*a*c**2*d*n*x*x**n/(
24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 3*a*c**2*d*x*x**n/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 36*a*c*
d**2*n**3*x*x**(2*n)/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 57*a*c*d**2*n**2*x*x**(2*n)/(24*n**4 + 50*n**3
+ 35*n**2 + 10*n + 1) + 24*a*c*d**2*n*x*x**(2*n)/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 3*a*c*d**2*x*x**(
2*n)/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 8*a*d**3*n**3*x*x**(3*n)/(24*n**4 + 50*n**3 + 35*n**2 + 10*n +
1) + 14*a*d**3*n**2*x*x**(3*n)/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 7*a*d**3*n*x*x**(3*n)/(24*n**4 + 50
*n**3 + 35*n**2 + 10*n + 1) + a*d**3*x*x**(3*n)/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 24*b*c**3*n**3*x*x*
*n/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 26*b*c**3*n**2*x*x**n/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) +
9*b*c**3*n*x*x**n/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + b*c**3*x*x**n/(24*n**4 + 50*n**3 + 35*n**2 + 10*
n + 1) + 36*b*c**2*d*n**3*x*x**(2*n)/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 57*b*c**2*d*n**2*x*x**(2*n)/(2
4*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 24*b*c**2*d*n*x*x**(2*n)/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 3
*b*c**2*d*x*x**(2*n)/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 24*b*c*d**2*n**3*x*x**(3*n)/(24*n**4 + 50*n**3
+ 35*n**2 + 10*n + 1) + 42*b*c*d**2*n**2*x*x**(3*n)/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 21*b*c*d**2*n*
x*x**(3*n)/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 3*b*c*d**2*x*x**(3*n)/(24*n**4 + 50*n**3 + 35*n**2 + 10*
n + 1) + 6*b*d**3*n**3*x*x**(4*n)/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + 11*b*d**3*n**2*x*x**(4*n)/(24*n**
4 + 50*n**3 + 35*n**2 + 10*n + 1) + 6*b*d**3*n*x*x**(4*n)/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1) + b*d**3*x*
x**(4*n)/(24*n**4 + 50*n**3 + 35*n**2 + 10*n + 1), True))

________________________________________________________________________________________

Giac [B]  time = 1.11209, size = 608, normalized size = 6.14 \begin{align*} \frac{24 \, a c^{3} n^{4} x + 6 \, b d^{3} n^{3} x x^{4 \, n} + 24 \, b c d^{2} n^{3} x x^{3 \, n} + 8 \, a d^{3} n^{3} x x^{3 \, n} + 36 \, b c^{2} d n^{3} x x^{2 \, n} + 36 \, a c d^{2} n^{3} x x^{2 \, n} + 24 \, b c^{3} n^{3} x x^{n} + 72 \, a c^{2} d n^{3} x x^{n} + 50 \, a c^{3} n^{3} x + 11 \, b d^{3} n^{2} x x^{4 \, n} + 42 \, b c d^{2} n^{2} x x^{3 \, n} + 14 \, a d^{3} n^{2} x x^{3 \, n} + 57 \, b c^{2} d n^{2} x x^{2 \, n} + 57 \, a c d^{2} n^{2} x x^{2 \, n} + 26 \, b c^{3} n^{2} x x^{n} + 78 \, a c^{2} d n^{2} x x^{n} + 35 \, a c^{3} n^{2} x + 6 \, b d^{3} n x x^{4 \, n} + 21 \, b c d^{2} n x x^{3 \, n} + 7 \, a d^{3} n x x^{3 \, n} + 24 \, b c^{2} d n x x^{2 \, n} + 24 \, a c d^{2} n x x^{2 \, n} + 9 \, b c^{3} n x x^{n} + 27 \, a c^{2} d n x x^{n} + 10 \, a c^{3} n x + b d^{3} x x^{4 \, n} + 3 \, b c d^{2} x x^{3 \, n} + a d^{3} x x^{3 \, n} + 3 \, b c^{2} d x x^{2 \, n} + 3 \, a c d^{2} x x^{2 \, n} + b c^{3} x x^{n} + 3 \, a c^{2} d x x^{n} + a c^{3} x}{24 \, n^{4} + 50 \, n^{3} + 35 \, n^{2} + 10 \, n + 1} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+b*x^n)*(c+d*x^n)^3,x, algorithm="giac")

[Out]

(24*a*c^3*n^4*x + 6*b*d^3*n^3*x*x^(4*n) + 24*b*c*d^2*n^3*x*x^(3*n) + 8*a*d^3*n^3*x*x^(3*n) + 36*b*c^2*d*n^3*x*
x^(2*n) + 36*a*c*d^2*n^3*x*x^(2*n) + 24*b*c^3*n^3*x*x^n + 72*a*c^2*d*n^3*x*x^n + 50*a*c^3*n^3*x + 11*b*d^3*n^2
*x*x^(4*n) + 42*b*c*d^2*n^2*x*x^(3*n) + 14*a*d^3*n^2*x*x^(3*n) + 57*b*c^2*d*n^2*x*x^(2*n) + 57*a*c*d^2*n^2*x*x
^(2*n) + 26*b*c^3*n^2*x*x^n + 78*a*c^2*d*n^2*x*x^n + 35*a*c^3*n^2*x + 6*b*d^3*n*x*x^(4*n) + 21*b*c*d^2*n*x*x^(
3*n) + 7*a*d^3*n*x*x^(3*n) + 24*b*c^2*d*n*x*x^(2*n) + 24*a*c*d^2*n*x*x^(2*n) + 9*b*c^3*n*x*x^n + 27*a*c^2*d*n*
x*x^n + 10*a*c^3*n*x + b*d^3*x*x^(4*n) + 3*b*c*d^2*x*x^(3*n) + a*d^3*x*x^(3*n) + 3*b*c^2*d*x*x^(2*n) + 3*a*c*d
^2*x*x^(2*n) + b*c^3*x*x^n + 3*a*c^2*d*x*x^n + a*c^3*x)/(24*n^4 + 50*n^3 + 35*n^2 + 10*n + 1)